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CDMS II data from the five-tower runs at the Soudan Underground Laboratory were reprocessed 
with an improved charge-pulse fitting algorithm. Two new analysis techniques to reject surface- 
event backgrounds were applied to the 612 kg days germanium-detector weakly interacting massive 
particle (WIMP)-search exposure. An extended analysis was also completed by decreasing the 
lOkeV analysis threshold to ~5keV, to increase sensitivity near a WIMP mass of 8GeV/c^. After 
unblinding, there were zero candidate events above a deposited energy of 10 keV and six events in the 
lower-threshold analysis. This yielded minimum WIMP-nucleon spin-independent scattering cross- 
section limits of 1.8x10“'*'* and 1.18x10“'** cm^ at 90% confidence for 60 and 8.6GeV/c^ WIMPs, 
respectively. This improves the previous GDMS II result by a factor of 2.4 (2.7) for 60 (8.6) GUV/c? 
WIMPs. 

PACS numbers: 95.35.+d, 95.30.Cq, 85.25.0j, 29.40.Wk 


I. INTRODUCTION 

The mass balance of the Universe is the subject of in¬ 
tense research and debate. Discrepancies between grav¬ 
itationally determined galaxy-cluster masses and their 
observed luminosities provided the earliest motivation 
for dark matter [1, 2]. Modern galactic rotation curves 
sharpen the argument [3-5], as do more recent studies of 
galaxy cluster dynamics [6, 7]. Likewise, spectroscopy of 
intergalactic x-ray-emitting gas [6, 7] and gravitational 
leasing [8-12] elevate these presumed mass discrepancies 
to the level of a crisis. 

To bring the data sets into consistency requires ei¬ 
ther modifications of gravity [13, 14] and/or large quan¬ 
tities of nonluminous matter. Observations of collid¬ 
ing clusters [15, 16] provide evidence for excess nonlu¬ 
minous matter, though alternate models of gravity and 
some nonluminous matter can apparently also reproduce 
such results [17, 18]. Nonrelativistic (cold) relic particle 
dark matter alone could resolve these discrepancies and 
is also considered an essential ingredient in gravitational 
simulations of the large-scale structure of the Universe. 
For example, the Via Lactea and Millennium simulations 
show excellent agreement with the observed large-scale 
structure of our Universe when cold dark matter is in¬ 
cluded [19-21]. 

While the observed galactic dynamics and large-scale 
structure naturally lead us to consider cold particle dark 
matter, cosmological measurements have been important 
in constructing a consistent model for the evolution of 
the Universe using such cold dark matter (CDM). The 
accelerating expansion of the Universe [22-24], big bang 
nucleosynthesis [25], baryon acoustic oscillations [26], 
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and the cosmic microwave background [27, 28] support 
a cosmology whose dominant components are dark en¬ 
ergy (which could correspond to a cosmological constant 
A) and nonbaryonic CDM. When interpreted within the 
framework of the ACDM model, these cosmological mea¬ 
surements enable precise determination of the CDM and 
dark-energy content of the Universe [29]. 

Particle physics provides clues as to the possible iden¬ 
tity of nonbaryonic CDM. It was realized early on that 
weakly interacting massive particles (WIMPs) with GeV- 
to TeV-scale masses could thermally freeze out in the 
early Universe to give the correct relic density [30]. Su¬ 
persymmetry provides a WIMP candidate in the light¬ 
est supersymmetric particle and has many other benefits, 
like solving the hierarchy problem [31]. Still other parti¬ 
cle physics models considered more recently [32] provide 
motivation for light WIMPs and linkages to the mat¬ 
ter/antimatter asymmetry of the Universe. 

WIMPs can be searched for directly through their scat¬ 
tering off nuclei in a terrestrial detector. Since these in¬ 
teractions are expected to be rare it is important for a 
dark-matter detector to have a low threshold (lOkeV or 
below) and excellent background rejection capabilities. 
The Cryogenic Dark Matter Search (CDMS) collabora¬ 
tion has developed cryogenic semiconductor detectors fo¬ 
cusing on those properties for the purpose of measuring 
WIMP-scattering events. 

This paper explores new analysis techniques using data 
taken by the CDMS experiment for the direct detection 
of WIMPs during the CDMS II running period [33-35] . 
There are four data periods associated with the data set 
used in this work, varying from 1 to 6 months in du¬ 
ration. The experiment used Ge and Si detectors and 
was located in the Soudan Underground Laboratory at a 
shielding depth of 2090 meters water equivalent (m.w.e.). 
The largest payload consisted of a total of 19 Ge (^240 g 
each) and 11 Si (~110g each) detectors. All detectors 
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are cylindrical, 7.6 cm in diameter, ^Icm in height, and 
arranged in five vertical stacks (“towers”) each includ¬ 
ing six detectors [36, 37] (see Fig. 1). The detectors are 
labeled TxZy where x (1-5) is the tower number and 
y (1-6) indicates the position within the stack (from top 
to bottom). 
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FIG. 1. (Color online) Isometric representation of the 
CDMS II detector arrangement and tower occupation, with 
direction of north indicated. The ZIP detectors in each tower 
are numbered 1-6 and are either silicon (yellow) or germa¬ 
nium (brown). The identification numbers of each detector 
were given to track the raw material the detectors were pro¬ 
duced from and the processing to which they were subjected. 
For example the detector “G9” is a germanium detector with 
identification number 9. Green stars indicate detectors that 
were used in this data analysis (see Sec. II). 

In each detector, particle interactions produce 
electron-hole pairs (ionization) together with phonons. 
The charge carriers are drifted by a small electric field 
(3V/cm for Ge and 4V/cm for Si) and collected on con¬ 
centric aluminum electrodes deposited on one of the flat 
faces of the crystal. On the opposite face, supercon¬ 
ducting transition-edge sensors (TESs) arranged in quad¬ 
rants collect phonons before thermalization. Charge and 
phonons are measured independently for the purpose of 
event-by-event discrimination; WIMP signal events will 
produce a nuclear recoil (NR) in the detector, whereas 
most background processes produce an electron recoil 
(ER). The ratio of ionization to phonon signal amplitudes 
allows discrimination of NRs from the far greater number 
of ERs with a rejection factor > 10^ [33] in the 10-100 keV 
region. Lead shielding further reduced gamma-induced 
backgrounds, and neutron-induced NR backgrounds were 


reduced with polyethylene shielding. The shielding is 
surrounded by an active scintillator veto to tag events 
induced by cosmic-ray muon showers. 

Recoils within ^10 fim of the detector surface can have 
reduced ionization signals because of poor charge col¬ 
lection, which can be sufficient to misclassify a surface 
ER as an NR. These reduced-charge ERs originate from 
three major sources: electrons produced by /3 emitters 
on or near the detector surfaces, electrons ejected from 
photons scattering in nearby material, and photons that 
interact near the detector surfaces. All of these mecha¬ 
nisms can lead to events in which an energy deposition 
is observed in only a single detector, producing “surface 
events” that can be mistaken for single-interaction NRs. 
The long-lived ^^^Rn daughter ^^°Pb, implanted in the 
detector surfaces and their copper housings, is the pri¬ 
mary source of such surface events. This class of events 
constitutes the dominant background for the CDMS II 
WIMP search and presents a considerable challenge [37]. 
Fortunately, recoils near detector surfaces are character¬ 
ized by prompt phonon absorption in the TESs. Con¬ 
sequently, the phonon signals for surface ERs are (on 
average) faster than for NRs in the detector’s bulk, en¬ 
abling surface-event background discrimination based on 
phonon-pulse timing. The event selection criteria derived 
from phonon timing (called “timing cuts”) presented here 
are essential to obtaining optimal WIMP sensitivity for 
the CDMS II data set because they mitigate the surface- 
event background, improving the overall ER rejection to 
> 10 ®. 

Nonsignal NR events have two known sources: Ge re¬ 
coils induced by scattering neutrons and nuclei 

from the decay of ^^°Po near the detector surfaces. The 
expected neutron background for CDMS II is roughly 
equal parts cosmic-ray muon-induced neutrons and radio¬ 
genic neutrons from trace contaminants in the shielding 
and detector. Radiogenic neutrons produced outside the 
shielding have too little energy to cause a detectable NR 
after penetrating the polyethylene. By combining simu¬ 
lations with in situ data, we estimate the neutron back¬ 
ground to be subdominant compared to the surface-event 
background (see Secs. VIIB and VII C). The ^°®Pb back¬ 
ground contribution is also subdominant in this analysis. 

The data set used in this work comes from the Ge de¬ 
tectors in the final CDMS II five-tower exposure and was 
acquired between July 2007 and September 2008. The 
total raw exposure for this running period was approxi¬ 
mately 612 kg days. The original analysis of this data set 
provided world-leading sensitivity to spin-independent 
elastic WIMP-nucleon scattering in 2010 when it was first 
published [33]. Here, we reevaluate this data set using im¬ 
proved data reduction algorithms and surface-event rejec¬ 
tion methods that reduce the expected surface-event con¬ 
tamination (“leakage”) in the WIMP signal region com¬ 
pared with the original analysis. In the 2010 analysis a 
lOkeV [38] analysis threshold was used to limit the ex¬ 
pected background to less than one event for the entire 
exposure. Intriguing results at low WIMP mass [39-43] 
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motivated us additionally to examine the data with re¬ 
duced thresholds. 


II. CDMS II DETECTOR PROPERTIES 


The CDMS II Z-sensitive ionization and phonon de¬ 
tectors (ZIPs) are operated at a temperature of ~50mK 
and feature six readout channels: two charge electrodes 
on one side and four phonon sensors on the opposite 
side [44] . The analysis undertaken here includes only Ge 
detectors, for which the ionization channels are biased 
with a 3 V potential across the crystal. Furthermore, five 
of the 19 Ge detectors were omitted from the analysis 
because of readout channel failures across all data sets. 
This gives a direct analogy to the original analysis of this 
data set, which used the same detector subset. Ioniza¬ 
tion channels are read out by a low-noise junction field- 
effect transistor circuit with an operating temperature of 
^150 K [45]. Phonons are detected by quasiparticle-trap- 
assisted electrothermal-feedback transition-edge sensors 
(QETs). The QET signal is amplified by supercon¬ 
ducting quantum-interference devices (SQUIDs) that 
are thermally coupled to the cryostat’s 600 mK cold 
stage [45]. Eigure 2 shows a schematic of the ZIP chan¬ 
nel layout. The outer charge electrode acts as a veto 
against events that deposit energy near the sidewalls, 
where higher background is expected and charge collec¬ 
tion is more likely to be incomplete. 


a) 

b) 



FIG. 2. (Color online) a) Inner and outer electrodes for mea¬ 
suring the ionization signal. The inner electrode extends from 
the center to a radius of 34.5 mm and the outer electrode ex¬ 
tends from a radius of 35.5 mm to just before the edge (0.5- 
2 mm; the diameter of the detectors is 76 mm). b) Phonon 
channels (QETs) A-D arranged in quadrants on the opposite 
side of the detector. Each phonon channel is composed of 
approximately 1000 TESs formed from thin films deposited 
onto the crystal surface and then photolithographically struc¬ 
tured. The TESs are connected in parallel to form the QETs 
and each QET is read out by its own circuit containing a 
SQUID array. 

Data taken during the CDMS II experiment are com¬ 
posed of calibration (using a ^^^Ba gamma source, or a 


^®^Cf neutron source) and WIMP search. The ^^^Ba cali¬ 
bration data were used to determine the energy scale and 
to study various systematic effects. The ^®^Cf calibration 
data were used to study detector response to NRs and act 
as a proxy for determining WIMP acceptance. 

The CDMS II front-end electronics issued experimental 
triggers in response to activity in either the surrounding 
scintillator veto or the ZIPs. For veto-triggered events, 
the front-end electronics required two or more veto pan¬ 
els to have coincident signals in excess of their hardware 
thresholds. ZIP-triggered events occurred when a com¬ 
posite phonon signal for any ZIP exceeded its preset dis¬ 
criminator threshold. The composite phonon signal was 
the analog sum of the four phonon signals from the ZIP 
with a 900-18000 Hz band-pass filter applied to reduce 
noise-like fluctuations. During ^^^Ba calibration runs, 
selective readout was employed because of the high rate; 
i.e. only detectors in the tower in which a ZIP trigger 
occurred were read out. During ^®^Cf calibration and 
WIMP-search runs, all detectors were read out in re¬ 
sponse to either veto or ZIP triggers. Veto-triggered 
WIMP-search events were recorded for use in studies of 
cosmogenically induced neutron backgrounds. Each trig¬ 
gered event includes two charge and four phonon traces 
(per ZIP readout). Each trace consists of 2048 digitized 
amplitudes acquired at a rate of 1.25 MHz, with 512 sam¬ 
ples prior to the trigger time. 

Charge carriers moving in the electric held of the de¬ 
tector generate phonons with a total energy proportional 
to the potential difference traversed (Neganov-Luke ef¬ 
fect [46, 47]). These “Luke phonons” are typically bal¬ 
listic and add to the phonon signal generated by the pri¬ 
mary recoil. The recoil energy was constructed from 
the total phonon energy Ep by subtracting a term pro¬ 
portional to the ionization-derived recoil energy Eq (or 
“ionization energy”) to account for the Neganov-Luke ef¬ 
fect: 


Er — Ep - ~Eq, ( 1 ) 

where e is the elementary charge, V is the absolute value 
of the operating potential of the detector, and e is the av¬ 
erage electron-hole pair creation energy. In Ge the value 
of e is approximately 3.0eV/pair for ERs [48]. The “ion¬ 
ization yield” (or yield) is dehned as the ratio of ion¬ 
ization energy to recoil energy (EqjEj.) and we require 
this quantity to be unity for ERs (see below). An NR will 
produce less ionization than an ER of equal recoil energy. 
This effect is well known [49] and provides the basis for 
the ZIP detector’s primary method of ER/NR discrim¬ 
ination. Eor recoil energies between 10 and 100 keV in 
Ge, the NR ionization yield varies between 0.2 and 0.3. 

Using the ^^^Ba calibration data, the charge and 
phonon amplitudes were calibrated so that the energies 
of known gamma lines are reproduced and the yield is 
unity for ERs. Specifically, the charge channels were cal¬ 
ibrated using the 356 keV line after a ^10% correction to 
the ionization amplitudes to account for small systematic 
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variations with interaction location (within the detector). 
We do not fully understand the origins of the systematic 
variation with interaction location, though it is empir¬ 
ically robust. One possibility is that the detectors are 
not neutralized uniformly by the infrared (940 nm) LED 
that is activated between data-taking periods to remove 
trapped charge. The relative phonon-channel calibration 
was performed by scaling the phonon-channel amplitude 
fractions-the amplitudes of individual channels divided 
by the sum of the amplitudes-to have equivalent distri¬ 
butions. Finally, the summed phonon energy was cali¬ 
brated by requiring the ionization yield to be unity on 
average for ER events in the 65-100 keV region. 


III. RAW DATA REDUCTION AND THE 
REPROCESSING 

Our analysis parameters are calculated from digitized 
charge and phonon pulses for each ZIP detector using 
a set of pulse-processing algorithms. These algorithms 
distill timing and amplitude information from the four 
phonon and two charge traces (per ZIP). Surface events 
are rejected using timing quantities derived from the 
raw data, such as the rise time r and the delay tdei of 
the largest of the four phonon pulses with respect to 
the faster ionization signal. The amplitude of the inner 
charge-channel pulse gives a measurement of Eq and the 
amplitude of the analog sum of the four phonon pulses 
gives a measurement of Ep after calibration. These en¬ 
ergy variables are used to construct the ionization yield 
{y) and recoil energy {Ej.) described in Eq. (1). 

The relative charge-channel amplitudes (charge “par¬ 
tition”), relative phonon-channel amplitudes (phonon 
“partition”), and relative phonon pulse timing provide 
information about an event’s position within the detec¬ 
tor. These additional parameters are also calculated us¬ 
ing the output of the pulse-processing algorithms and are 
used for the phonon event-position-based correction and 
charge-derived fiducial-volume restrictions. 

This section gives an overview of the most important 
algorithms for constructing these analysis parameters. It 
also explains the upgrade to the charge-pulse processing 
that motivated the reprocessing of this data set. 


A. Parameter extraction 

Timing parameters. Timing estimators are derived 
using an algorithm that steps along a low-pass-filtered 
trace to identify the pulse rise time (RT) and fall time 
(FT). We call this the RT-FT-walk algorithm. The traces 
are first filtered with a low-pass Butterworth hlter [50]. 
Two sets of parameters are produced; one where the 
cutoff frequency is 50 kHz and one where it depends on 
the signal-to-noise of the trace. The filter removes high- 
frequency noise and effectively smooths the pulse for an 
improved determination of the RT and FT at various 


percentages of the pulse maximum. RT and FT informa¬ 
tion is determined by “walking” along the Hltered trace 
starting at the maximum and identifying the times at 
which the respective threshold levels are reached. For 
example, it infers the time at which the rising pulse edge 
reaches 20% of the pulse’s maximum amplitude. The RT- 
FT-walk algorithm was applied to all charge and phonon 
traces. The rise time r is computed as the 10-40% time 
span along the rising edge of the largest of a detector’s 
four phonon pulses using the RT-FT-walk algorithm. 

Optimal filtering. A pulse-template optimal hlter 
(OF) [51] is used to produce the best resolved energy 
quantities. The OF has superior energy resolution com¬ 
pared to pulse-integral quantities and is well suited to the 
analysis of small recoil energies [52]. A template for the 
expected pulse shape is ht to the pulse in Fourier space, 
deweighting the frequency bins with high noise. The 
frequency deweighting is done for each individual data 
“series” - a data-taking block normally lasting between 
10 and 12h-using noise power spectral densities (PSDs). 
The PSD is constructed from randomly triggered traces 
taken before the series to sample the noise environ¬ 
ment. Phonon pulse templates are two-exponential func¬ 
tional forms, with rise- and fall-time parameters tuned to 
match individual detectors by htting to an average pulse. 
Charge pulse templates were produced empirically by av¬ 
eraging normalized data traces. This htting procedure is 
done for all possible pulse template delays and the best- 
ht delay is chosen. 

A single-pulse-template OF is performed on each 
phonon pulse separately and on the sum of a detector’s 
four phonon traces. For charge pulses, to account for 
crosstalk between the inner and outer channels, the two 
pulses are ht simultaneously with a crosstalk-correcting 
OF (OFX) described in more detail below. During the 
OF ht an array of best-ht amplitudes is produced that 
corresponds to each possible delay of the template with 
respect to the experimental pulse. The maximum such 
amplitude in a preselected delay window is chosen and 
its time dehnes the delay quantity. For an OF using 
a single pulse template this maximum amplitude choice 
also has the lowest in the preselected window. The 
OF amplitude for the summed phonon trace gives Ep af¬ 
ter calibration and the inner charge amplitude from the 
OFX procedure gives Eq after calibration. The charge- 
to-phonon delay tdei is determined from the difference 
between the OFX charge delay and the 20% crossing of 
the largest phonon signal as determined by the RT-FT- 
walk algorithm. The delay is computed as follows: 


^del — ^p20 ~ ^OFX, (2) 


with toFX being the OFX ionization delay relative to the 
global trigger time and tp 2 o being the 20% crossing time 
of the largest phonon signal relative to the global trigger. 
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B. Parameter corrections 

Phonon position correction. The phonon sensor re¬ 
sponse is highly position dependent [53] . In order to opti¬ 
mize the derived phonon variables, a position correction 
was used to modify energy, timing, and yield quantities 
by normalizing to the mean of nearby events. We first 
defined averaging neighborhoods using a five-dimensional 
metric that included total phonon energy, two relative 
phonon delay parameters, and two phonon energy par¬ 
tition parameters [54]. For each parameter to be cor¬ 
rected, a lookup table containing the neighborhood aver¬ 
ages for each five-dimensional bin was then constructed 
using bulk ER events confirmed to be of good quality 
from ^^^Ba calibration. To apply the correction we cre¬ 
ated a corrected version (vc) of a given parameter (v) as 
in Eq. (3): 

Vc = -yv (■y) global, (3) 

where the subscript “global” refers to averaging over the 
whole calibration data set, and “bin” refers to averag¬ 
ing over one neighborhood. After the lookup table was 
applied to the ^^^Ba calibration data, the selection of 
bulk ERs was refined to create the final correction ta¬ 
ble, which was then applied to the WIMP-search data. 
The final energy, timing and yield quantities were made 
more uniform across the detector for bulk events by us¬ 
ing this procedure, thus improving the rejection of surface 
events. 

Charge crosstalk. Charge signal quantities use the 
OFX algorithm to account for crosstalk. Using separate 
pulse-template OFs for the inner and outer electrodes, 
a single delay was chosen. In the original analysis of 
our data set [33], this delay corresponded to the delay 
that made the sum of the two OF amplitudes maximal. 
While this is likely to be very close to the delay that 
minimizes the of th® fif, it is not guaranteed because 
of the combined fitting of the inner and outer electrodes. 

C. Reanalysis motivation 

After finalizing the first analysis of this data set in 
2009 [33] we noticed that one of the WIMP candidate 
events had an OFX delay value that did not correspond 
to a global minimum in the of the fit (see Fig. 3). This 
issue was one of the main reasons for proceeding with a 
reanalysis. 

The reason for implementing the maximum-amplitude 
method instead of a full x^ minimization in the origi¬ 
nal analysis was to save computing time; a global trigger 
issued by one detector forces the readout of all detec¬ 
tors, giving a substantial number of traces with near-zero 
amplitudes that need not be processed with the time- 
consuming minimization. The result of using this 
maximum-amplitude method is to smear the ionization- 
to-phonon delay distributions slightly, which can cause 



FIG. 3. (Color online) The OFX as a function of delay 
for the T3Z4 WIMP candidate from the original analysis [33]. 
The delay chosen by the original ionization fit to this event 
(red dashed line) is not a global minimum in y^, rather it is 
displaced ~4 /rs from the global minimum (green dot-dashed 
line). While this small error in the inferred delay for this and 
other events did not affect the measured energy significantly 
(a ~2% decrease), it caused an error in the surface-event- 
rejection timing parameter for low-energy events. 

additional background events to leak past the timing 
cuts. This delay smearing affects lower-ionization events 
more than higher-ionization events. In particular, NRs 
that deposit energy in a single detector (“single scat¬ 
ters”) are more susceptible than the (higher-ionization) 
low-yield ERs used to quantify the expected leakage of 
surface events into the signal region. 

The OFX inaccuracy made only a small contribution to 
the previously published limit, and background estimates 
were adjusted upward in the original result to account 
for the resulting larger uncertainties at low energies [33]. 
Nevertheless, the number and character of the candidate 
events provides useful information on the possibility of a 
signal. It is very likely that the candidate event featured 
in Fig. 3 would have been removed by the timing cut in 
the original work had the ionization delay at the global 
minimum of the y^ been chosen, but using the improved 
OFX procedure also has the potential to cause previously 
excluded events to become WIMP candidates. 

For the reasons stated above, the data were repro¬ 
cessed, selecting the global minimum of the y^ in the 
OFX algorithm for most ionization pulses, rather than 
the maximum-amplitude method. Because of the large 
number of traces in the raw data consistent with noise, 
traces with pulses corresponding to charge energies be¬ 
low a detector-dependent threshold (0.94keVee [38] on 
average across the 14 detectors used in this work) were 
still processed with the maximum-amplitude method in 
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order to save processing time. This has no effect on the 
results presented here because this energy is below our 
lowest analysis threshold. 


IV. DATA SELECTION AND EFFICIENCIES 

Once the data set was calibrated and position cor¬ 
rected, we created several selection criteria (or “cuts”) to 
produce the cleanest signal-region sample possible. Since 
the rejection of nonsignal events is typically not perfect, 
a signal event retention efficiency was computed for each 
cut. The cuts and efficiencies are covered in this sec¬ 
tion, starting with the efficiency for the trigger-a selec¬ 
tion criterion that is made in hardware before the events 
are recorded. The result is a signal-region sample corre¬ 
sponding to a well-known exposure. 

A trigger efficiency for each detector was determined as 
a function of phonon energy and converted to a function 
of recoil energy using our measured NR ionization yield 
from ^®^Cf data. A given detector’s trigger efficiency 
was calculated using all WIMP-search events in which 
another detector caused an experimental trigger. When 
the global trigger initiates a readout of all the detectors, 
other delayed instances where phonon pulses were in ex¬ 
cess of the corresponding discriminator thresholds were 
recorded into a logical buffer. Regardless of the con¬ 
tent of these trigger records, for each event and detector 
the optimal-filtering techniques described in Sec. II were 
used to reconstruct the total phonon energy (Ep) from 
Eq. (1). The trigger efficiency curve is then the ratio of 
two recoil-energy spectra: the spectrum of events with a 
phonon trigger in both the detector in question and an¬ 
other detector divided by the spectrum of events with a 
phonon trigger in another detector. As shown in Fig. 4- 
combined across the 14 Ge detectors considered here-the 
trigger efficiency is ~100% for recoil energies above our 
standard analysis threshold of 10 keV. Note that the ex¬ 
tended analysis described in Sec. V D requires the trigger 
efficiency down to ~5 keV. 

Events above the trigger threshold were subjected to 
five classes of cuts used to isolate high-quality signal can¬ 
didates in the WIMP-search data: cuts that remove time 
periods of reduced overall data quality; event-level data- 
quality cuts; an ionization-based fiducial-volume cut to 
reject events near the detector sidewall; an ionization- 
yield cut; and one of three phonon timing cuts. Cuts 
that remove time periods cause a loss of experimental live 
time because events are removed uniformly for all detec¬ 
tors and in all kinematic variables. After the live-time 
cut, the total remaining live time for each detector was 
computed and stored. All other cuts cause reductions in 
NR detection efficiency, some of which vary significantly 
as a function of recoil energy. We computed these “effi¬ 
ciency functions” and applied them to each detector in 
conjunction with the live time to compute the final ex¬ 
posures. 

The NR detection efficiency functions were estimated 


C 

a> 

o 

0) 


- trigger 


ri-.-rT 

1-1-ri- 

n ■ 

TT. 

r-rn-r 

n-r-r 

p 








= 

: 

uati 

ty 

ldu< 

>ial ' 

fOltl 




i 

: i/” ... 

: i /’ . 


fn© 

. 


i 

1. I'' 




V 

ield 



■ -I 









i 

^ 1 / 




^ 




i 

- f 1 / 





1 

>cl-5C 



: i|' 

■ i 








i 

C J i 
c if! 








i 



1 ■ ■ I ■ 

. 

. 

. 

. 

. 



recoil energy [keV] 


FIG. 4. (Color online) Detection efficiency combined across 
all detectors as a function of energy for different classes of 
cuts. The curves show the total efficiency as more cuts are 
added beyond the hardware trigger (pink triple-dot-dashed): 
event-level data-quality cuts (green solid), ionization-based 
fiducial-volume cut (magenta dotted), ionization-yield cut 
(maroon dot-dashed), and 5d-x^ (see Sec. V C) phonon tim¬ 
ing cut (blue dashed). The falloff of the data-quality cut 
efficiency below ~20keV is due to the interplay between the 
ionization threshold and the requirement that the ionization 
yield be Scr below the ER band (see main text). The vertical 
black dashed fine is our standard lOkeV analysis threshold. 

in bins of recoil energy. Well-motivated functional forms 
were fit to the bin-wise estimates, with the best-fit re¬ 
sults shown in Fig. 4. Most cuts have little effect on 
the acceptance of nuclear recoils. The ionization-based 
fiducial-volume cut and the phonon timing cuts cause 
the greatest loss of signal acceptance. Of the three tim¬ 
ing cuts described in Sec. V, use of the “5d-x^” timing 
cut results in the highest final acceptance above lOkeV, 
about 50% at ~35 keV. The 5d-x^ timing-cut efficiency is 
shown in Fig. 4 as an example while the other timing-cut 
efficiencies are detailed in Sec. VE. 


A. Live-time cuts 


Live time was removed during periods with disabled 
readout channels, poor detector neutralization (charac¬ 
terized by increased levels of charge trapping in the crys¬ 
tal bulk [52]), decreased resolution, improper experimen¬ 
tal configurations, and trigger anomalies that consist 
of isolated bursts of events or incorrect phonon trigger 
threshold. The Soudan Underground Laboratory also 
houses a neutrino detector, MINOS, to measure proper¬ 
ties of the “NuMI” neutrino beam originating from Fermi 
National Accelerator Laboratory [55]. While it is very 
unlikely for CDMS to observe any beam-induced events. 






































the loss of live time incurred by removing all time pe¬ 
riods coincident with a NuMI beam spill is negligible. 
The NuMI neutrino-beam cut is implemented as a re¬ 
striction on the live time, removing events within 60 /iS 
of the arrival of the neutrino beam’s 10 /iS spills. After 
application of these cuts, 612 kg days of total exposure 
remain. We define this as our “raw exposure.” Out of 
a total of 36% loss of live time (64% data-taking efh- 
ciency) the largest contributors were bad environmen¬ 
tal configurations (~15%), failed KS tests (~8%), cry- 
ocooler noise (^6%), poor phonon reconstruction (~4%), 
and fundamental hardware failures (~3%). Bad environ¬ 
mental configurations can consist of many effects includ¬ 
ing trigger bursts, high charge noise, or insufficient LED 
flashing. 


B. Quality cuts 

The event-level class of quality cuts consists of sev¬ 
eral components. A “glitch” cut removes events that 
have phonon pulses resembling electronic noise in the 
phonon readout chain. These events are typically cor¬ 
related across multiple detectors and are characterized 
by phonon pulses with fall times shorter than the time 
scale expected for phonon dissipation. Such events are 
less likely to appear in the ionization readout chain; con¬ 
sequently, they are effectively removed with a cut that 
rejects events in which the phonon pulse multiplicity is 
significantly larger than the ionization pulse multiplic¬ 
ity [56] . The muon-veto cut uses the 2 in. thick scintil¬ 
lator panels that surround the CDMS II experiment and 
the following two rejection criteria: an event was removed 
if 1) a ZIP event is associated with a muon-like energy 
deposition (0.58 V or ~3.8MeV deposited) in any veto 
panel occurring between 185 /iS before and 20 fis after the 
ZIP trigger time, or 2) there is any veto activity above 
the scintillator-panel hardware threshold (~0.23V) and 
within the 50 /iS before a ZIP trigger. Both criteria are 
used for the muon-veto cut because it makes the cut 
stronger and has only about a 2% efficiency loss across 
all recoil energies. Finally, because WIMPs will not in¬ 
teract in more than one detector within a given event, we 
defined a “single-scatter” cut to select events involving a 
single ZIP detector as follows: 1) the phonon signal must 
be greater than six standard deviations above the mean 
electronic-noise level in the detector under consideration; 
and 2) the signal in all other detectors must be within 
four standard deviations of the means of their respec¬ 
tive noise distributions. Ionization channels were used 
for multiple-scatter rejection for those detectors with de¬ 
graded phonon channel performance. Thus, Ge and Si 
detectors that are not part of the WIMP-search expo¬ 
sure were still live with respect to identification and re¬ 
jection of multiple-scatter events. The single-scatter cut 
efficiency for each detector was estimated from the frac¬ 
tion of randomly triggered events (i.e., electronic noise) 
that satisfy criterion 2). The combined efficiency of the 


glitch, muon-veto and singles cuts varies by detector and 
over time, and ranges from 96% to 98%. Another aspect 
of the singles cut is its use to select certain samples from 
calibration data. For example, multiple-scatter events 
that have their secondary scatter in the detector above or 
below the triggering detector are said to be “face tagged.” 
For these events the primary recoil is biased toward the 
direction facing the multiples tag. When using the mul¬ 
tiples in background estimations, separating the events 
by face helps decrease the systematic uncertainties (see 
Sec. VII A). Sometimes events with multiples toward the 
phonon readout side are called “phonon-side” and those 
with multiples toward the charge readout side are called 
“charge-side.” 

Some quality cuts depend explicitly on the kinematic 
variables. One of these is an ionization threshold cut, 
which requires events to have reconstructed ionization 
signals greater than 4.5 standard deviations above the 
mean noise level. This cut removes nearly all events 
with zero charge collected (those from very near the side 
walls), resulting in <0.1 of such events over the whole 
exposure. A charge-pulse reconstruction-quality cut was 
also defined, requiring the value of the OFX fit to be 
less than an energy-dependent threshold. This cut se¬ 
lects events with high-quality charge energy estimators, 
suppressing those with excess electronic noise or the oc¬ 
currence of multiple pulses within a single event trace 
(referred to as “pile-up”). Its efficiency was measured as 
a function of charge energy using ^^^Ba calibration data 
and translated into recoil energy using the average NR 
ionization yield measured from ^®^Cf calibration. Above 
~20keVee [38] the efficiency is constant and > 99%. It 
decreases slightly at lower recoil energies. 

Charge collection is reduced near the cylindrical walls 
of the detector, causing lower ionization yield and thus 
poor NR to ER discrimination. The ionization-based 
fiducial-volume cut accepts only events with an outer 
electrode signal consistent with noise. The efficiency of 
this cut (see Fig. 4) was estimated using the ^®^Cf cal¬ 
ibration data, including a small correction (^13% max¬ 
imum across the energy range) accounting for residual 
leakage of ERs into the NR signal region in these data. 
Another ^5% correction was applied across the whole en¬ 
ergy range to account for detector self-shielding to neu¬ 
trons and multiple scattering; this correction was based 
on a Monte Carlo simulation of neutron scattering. 


C. Yield cut 

An energy-dependent cut on ionization yield is used as 
the primary method for discriminating NRs from ERs. 
An NR “band” was derived for each detector by fitting 
the distribution of NR yields in ^®^Cf calibration with a 
Gaussian hypothesis in bins of recoil energy. The collec¬ 
tions of best-fit Gaussian means and widths were then fit 
with energy-dependent functional forms, giving smooth 
parametrizations versus recoil energy of the average yield 
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and the yield resolution for NRs. The functional form 
for the means was inspired by the Lindhard theory [57] 
{y = a-E^, where a and b are fitted parameters), whereas 
the widths are fit to a power law below a fitted energy 
threshold and a constant above [58] . ER bands were sim¬ 
ilarly constructed from ^^^Ba calibration. The primary 
ionization-yield cut requires that events be located within 
the ±2cr width of the NR band. By construction, the 
selection efficiency is ~95% and roughly constant with 
energy (see Fig. 4). Variations in detector response over 
the course of the WIMP search caused slight variations 
in each detector’s NR and ER bands. Consequently, the 
bands depend on both detector and time. The time vari¬ 
ation is represented by widths of the band lines shown in 
Fig. 5, where they are plotted with the ^®^Cf data from 
which the NR bands are derived. 

To prevent ERs from entering the signal region at low 
recoil energy where the ER and NR bands overlap (see 
Fig. 5), the yield-based discrimination is refined to in¬ 
clude a requirement that candidate events lie at least 3cr 
below the mean of the ER band (see Sec. IV D). This 
condition and the ionization threshold have the greatest 
impact on the overall NR detection efficiency at low en¬ 
ergies. The efficiencies for these two cuts are 100% for 
recoil energies greater than ^15-25 keV (depending on 
the detector) and decrease rapidly to zero for lower en¬ 
ergies. This behavior is illustrated in Fig. 4, where the 
efficiency of these two cuts is shown combined with the 
other event-level quality cuts. 



FIG. 5. (Color online) ER and NR bands defined in the 
plane of yield versus recoil energy for a representative detector 
(T1Z2). The ER band mean ±3 (j curves (red hatched), NR 
band mean ±2(j curves (blue fine cross-hatched), ionization 
threshold 3.On and 4.5(7 curves (green coarse cross-hatched), 
and the 10 keV threshold line (black dashed) are shown super¬ 
imposed onto ^®^Cf data (points). The widths of the bands 
represent the time-period variation of the fits that are used 
to define them. 


D. Data regions 

To assess the effectiveness of the timing cuts defined 
in Sec. V, several standard data regions in the nontim¬ 
ing variables were defined. This terminology will be used 
throughout the rest of this work to describe data sam¬ 
ples used for timing-cut tuning and sideband-style back¬ 
ground estimations. All of the event samples include the 
data-quality and fiducial-volume restrictions. Further, 
they include the requirement that events have charge 
energies greater than 4.5 standard deviations above the 
mean electronic-noise level. 

1. Nuclear-recoil single scatter (NRSS). Events 
that are below the ER band mean -3cr line, within 
the NR band mean ±2(7 lines, and are single scat¬ 
ters. The portion of the NRSS events that pass the 
timing cut and are in the range of 10-100 keV make 
up the WIMP signal region. The Sd-y^ timing cut 
uses a slightly modified signal region in that the 
width of the NR band is optimized. We use the 
term NRSS there as well, leaving the precise width 
of the NR band to be determined from context. 

2. Nuclear-recoil multiple scatter (NRMS). 

Events that are below the ER band mean -3tT line, 
within the NR band mean ±2cr lines, and are not 
single scatters. 

3. Wide-band (WB). Events that are below the ER 
band mean -5a line, and above the NR band mean 
±2(7 line. 

4. Wide-band multiple scatter (WBMS). Events 
that are below the ER band mean -5a line, above 
the NR band mean ±2(7 line, and are not single 
scatters. These events are a good representation of 
surface events. 


V. PHONON TIMING DISCRIMINATION 

ZIP detectors have an excellent ability to discriminate 
NRs from ERs if the energy depositions occur away from 
the surfaces of the detector. But background surface 
events can populate the NR band, despite being ER in 
nature, since they can have reduced ionization yield. We 
remove surface events from our WIMP candidate sample 
using the timing characteristics of the phonon signals. A 
small number of surface events, however, can survive into 
the signal region. We call these leakage events. The def¬ 
inition of the timing cut affects the surface-event leakage 
and the total exposure (through the NR acceptance), so a 
crucial element in the timing-cut construction is tuning 
for optimal WIMP-detection sensitivity. Three timing- 
cut constructions are reviewed in this section. 

Having three independent methods for surface-event 
rejection gives a handle on the systematic uncertainty 
of the leakage estimates. The WIMP limits for these 
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three realizations of surface-event rejection and sensi¬ 
tivity maximization are presented in Sec. VIII. Table I 
gives a summary of the different timing parameters used 
for each timing cut construction: the “classic,” “neural- 
network,” and “5d-x^” analyses. The choice of param¬ 
eters used in particular timing-cut constructions is ex¬ 
plained in the following sections. 


Quantity 

Description 

Analysis 

r 

VF phonon rise time 

NN, 5d-x" 

^del 

VF phonon delay 

NN, 5d-x^ 

r 

CF phonon rise time 

Classic, 5d-x^ 

^del 

CF phonon delay 

Classic, 5d-x^ 

■froTO 

CF phonon 40-70% rise time 

5d-x' 

w 

CF phonon pulse width 

NN 

P5070 

phonon 50-70 kHz power 

NN 


TABLE I. Brief description of each of the phonon timing 
quantities used in this work. The abbreviation “VF” indi¬ 
cates the use of variable-frequency filtering prior to the appli¬ 
cation of the RT-FT-walk algorithm. “CF” indicates a con¬ 
stant (50 kHz) filter prior to that algorithm (see Sec. Ill A). 
Under “Analysis,” “NN” refers to the neural-network timing 
analysis. 


Each of the timing cuts was optimized to produce the 
best expected sensitivity to WIMPs given the expected 
leakage. Since we do not know the WIMP mass, a “tar¬ 
get” value is chosen for each analysis (60 GeV jc? for 
most analyses in this work) and the expected sensitiv¬ 
ity is maximized given that WIMP mass. The spectrum- 
averaged exposure (SAE) is a way to quantify the amount 
of the raw exposure (MT = detector Mass x live Time) 
that is utilized toward the WIMP search over the analy¬ 
sis energy range, given a WIMP recoil spectrum for mass 
of f{Er;m^). The SAE is computed as follows: 


El, Eh) =MT 


/I'- dE,e{E,)f{E,;m^) 

Sei dErf{Er-,m^) 


(4) 


where Er is the recoil energy, Ei {Eh) is the lower (up¬ 
per) signal-region energy limit, and e{Er) is the cumu¬ 
lative signal acceptance efficiency for all cuts at a recoil 
energy Er- Note that the SAE is equal to the raw expo¬ 
sure only when the analysis efficiency is unity over the 
entire energy range and that only SAEs with the same 
WIMP-mass assumptions and signal-region energy range 
are strictly comparable. Sometimes the right-hand side 
of Eq. (4) with MT divided out is called the spectrum- 
averaged efficiency. The SAE is computed on a detector- 
by-detector basis and then summed in the final analysis. 

Before looking at the events in the final signal region 
(see Sec. VI on “unblinding”) we used the expected leak¬ 
age and the SAE of each timing cut optimization to calcu¬ 
late the expected sensitivity. The expected sensitivity is 
computed by using the expected leakage and calculating 
the upper limit of counts at the 90% C.L. This is normal¬ 
ized by the SAE to produce the lowest WIMP rate the 
experiment is sensitive to. For the 10 keV threshold anal¬ 
ysis the 5d-x^ timing cut had the best sensitivity, while 


our lower-threshold analysis showed the classic method to 
have the best sensitivity. The main results of this work 
are therefore the limits of the 5d-x^ analysis for WIMP 
masses above I1.3GeV/c^, and the limits of the classic 
method for WIMP masses below lI.3GeV/c^. 


A. Classic timing analysis 


The phonon timing cut strategy that was used in 
the original analysis of these, as well as earlier CDMS 
data [33, 36, 59] was also used in our reanalysis and pro¬ 
vides a point of comparison between the two. 

Our “classic” timing parameter is defined as the sum of 
two quantities: the delay (tdei) and the 10-40% rise time 
(r), both derived from the most energetic phonon signal 
among the four sensors (see Sec. Ill) . The sum is approx¬ 
imately the optimal combination of these two variables, 
as can be seen in Fig. 6. A timing cut is defined as a set 
of detector-dependent thresholds on the distributions of 
this parameter, below which all events are rejected. The 
thresholds were determined by an optimization scheme 
that approximately maximizes the sensitivity to a WIMP 
with a mass of 60GeV/c^. Practically, this was accom¬ 
plished by maximizing the WIMP-search exposure (as 
measured with ^®^Gf NRs) while keeping the total leak¬ 
age approximately equal to a “target” leakage of ~0.5 
events. This approximately maximized the 60 GeV jc? 
WIMP sensitivity while keeping the total expected leak¬ 
age well under one event [60]. 

The expected surface-event leakage was estimated from 
representative ^^^Ba calibration and sidebands in the 
WIMP-search data that are insensitive to WIMPs. The 
surface-event background estimates for the three timing- 
cut strategies are described in detail in Sec. VII A. Al¬ 
though the quality cuts remove most unusual events, 
data-reconstruction artifacts occasionally result in events 
with extreme kinematic quantities (“outliers”). To pre¬ 
vent such outliers from skewing the timing-cut optimiza¬ 
tion, a consistency cut rejects events for which r -|- tde\ is 
greater than 32 ps or r-tdei falls outside the 0.5% and 
99.5% quantiles of the combined ^^^Ba and ^®^Gf data 
sets. Figure 6 shows the classic timing cut in the de¬ 
lay versus rise-time plane, and Fig. 7 shows the cut in 
the yield versus timing-parameter plane for an example 
detector [61]. Applying this timing cut results in a to¬ 
tal SAE summed over detectors of 220 kg days between 
10 and 100 keV for a 60GeV/c^ WIMP and an expected 
surface-event leakage of 0.64 (-1-0.17-0. 15)stat events, cal¬ 
culated after unblinding. The post-unblinding calcula¬ 
tion is generally more accurate because it makes use of 
the previously sequestered (see Sec. VI) NR single scat¬ 
ters that failed the timing cut (see Sec. VII). 
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FIG. 6. (Color online) Classic timing cut in the delay versus 
rise-time plane for a representative detector (T1Z5). ^®^Cf 
NRSS and NRMS events (blue open circles) and WBMS 
events (red crosses) are shown. Accepted events lie within the 
“consistency” region (black dashed lines) and to the upper 
right of the discrimination cut (black solid line). Considering 
the ^^^Ba WBMS events (red crosses), 2381 events fail the 
timing cut and four events pass. 
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timing parameter [|ds] 


FIG. 7. (Color online) Classic timing cut in the ionization- 
yield versus timing-parameter plane for a representative de¬ 
tector (T1Z5). ^®^Cf NRSS and NRMS events (blue open cir¬ 
cles), ^^^Ba WBMS events (red crosses), and ^^^Ba events in 
the ER ±2cr band (black filled circles) are shown. Accepted 
events lie to the right of the timing parameter line (black 
dashed line) and in the NR band (black solid lines). Consid¬ 
ering the ^®^Ba WBMS events (red crosses), 2381 events fail 
the timing cut and four events pass. 


B. Neural-network timing analysis 

A neural-network technique was used to develop a tim¬ 
ing cut using four timing parameters: the previously 
defined phonon delay fdei (i) and rise time r (ii); the 
phonon pulse width w (iii), defined as the time difference 
between the 80% points on the rising and falling edges 
of the largest-amplitude phonon pulse; and the spectral 
power of the largest phonon pulse P 5070 (iv), integrated 
between 50 and 70 kHz [62] . These parameters were cho¬ 
sen because they showed the most promising discrimina¬ 
tion in their one-dimensional distributions. 

These four variables were fed into a principal compo¬ 
nent analysis [63] . Principal component analysis is a sta¬ 
tistical method for determining a unitary transformation 
that takes N possibly correlated input vectors and re¬ 
turns N output vectors that are linear combinations of 
the input vectors (Af = 4 in this case). The output vec¬ 
tors are ordered by their statistical variance, so that the 
output vector with the ith-highest variance is called the 
ith principal component. Since the input vectors can 
have different characteristic scales and are dimensional, 
the input vectors were normalized to zero mean and unity 
variance so that the ordering of statistical variances of the 
output vectors is meaningful in an absolute sense. 

Neural-network computational complexity scales 
poorly with the number of input parameters. Therefore, 
only the first two principal components (i.e., those 
with the highest variance) were selected as inputs for 
the neural network. Given input parameters with 
similar intrinsic resolution and physical relevance, the 
high-variance combinations will be those that maximally 
separate the distinct populations (i.e., NRs and ER 
surface events). Principal components were selected 
separately for each detector and for several bins of 
time spanning the WIMP-search data set. The latter is 
necessary to capture changes in detector performance 
caused by variations in operating conditions. The 
principal component rotation showed that all four 
timing parameters contribute significantly to the first 
and second principal components in most cases. This 
generally indicates that the use of the extra parameters 
(as compared to the classic analysis) is beneficial even 
though in the end the sensitivity change is not dramatic 
(see Sec. VIII A). 

The neural network that was used is a multilayer per- 
ceptron with one hidden layer, 30 neurons, and a logistic 
sigmoid activation function. The NETLAB package [64] 
for MATLAB was used to perform this analysis. Train¬ 
ing samples for surface events and NRs were selected 
from the ^^^Ba and ^®^Cf calibration data, respectively. 
Events from the NR training sample were assigned a tar¬ 
get output value of 1, while surface events were assigned 
a target value of 0. Data were sorted into two bins of 
recoil energy, above and below 30keV, in order to take 
the energy dependence of the timing parameters into ac¬ 
count. Finer energy binning was not possible because of 
small statistics in the training samples. 
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Separate neural networks were trained for each com¬ 
bination of detector, time period, and energy bin us¬ 
ing the standard back-propagation-of-errors training al¬ 
gorithm [64]. Once the neural networks were trained, 
they assigned a numerical value to each event in the 
range 0-1. An output value close to 1 (0) corresponds 
to events of NR-like (surface-event-like) character. The 
distributions of this output parameter for the calibration 
data, as seen in Fig. 8, were then used to set a threshold 
for each detector, time period, and energy bin. 

The thresholds were set such that the WIMP-search 
exposure is maximized using a target surface-event leak¬ 
age of 0.5 events. As in the classic analysis, this 
target leakage approximately maximized the 60 GeV j(? 
WIMP sensitivity and the same optimization procedure 
was used. Following optimization, the total SAE for a 
60 GeV/c^ WIMP is about 216 kg days, with an expected 
leakage of 0.87 (-1-0.24-0.21)stat events. Similar to the 
leakage estimate for the classic timing cut, this expected 
leakage was estimated after unblinding for better accu¬ 
racy (see Sec. VII). 


C. 5d-y^ timing analysis 

The 5d-x^ surface-event rejection [65] was imple¬ 
mented by differentiating events based on a goodness of 
fit to two event-type hypotheses [66, 67]: Xn 
Xb for surface events. Five timing quantities were used 
to form each x^ value. Three of the quantities are mea¬ 
sures of rise time for the largest-amplitude phonon chan¬ 
nel: two based on the 10-40% rise time (r and r); and 
one based on the 40-70% rise time ( 74070 ). The remain¬ 
ing two quantities correspond to the delay of the phonon 
pulse relative to the prompt charge pulse - one computed 
using a variable and one a constant-frequency RT-FT 
walk (t(jei and tdei)- The inclusion of different measures 
of the same physical quantities (delay and rise time here) 
increases the robustness of the x^ value and is more ef¬ 
fective at identifying outliers. Therefore, the delay and 
rise-time parameters with good one-dimensional discrim¬ 
ination and the least redundancy (correlation) were cho¬ 
sen. 

Event samples of neutrons and surface events taken 
from calibration data were used to constrain the timing- 
quantity distributions for each event type. For each 
detector, calibration data were separated into neutron, 
charge-side surface-event, and phonon-side surface-event 
samples. The energy-dependent means, fi{Er;a)-a vec¬ 
tor of the timing-quantity distribution means for each 
event type a-were then fit to the empirically motivated 
functional form 

n{Er; a) = ai{a) -f a 2 {a)E^ + (5) 

where the ai(a) are free parameters for each timing quan¬ 
tity (the vector indices) and for particle type a. The 
term is observed to improve the fit. The covariance ma- 




transformed neural-network output 

FIG. 8. (Color online) Neural-network training distributions 
of a representative detector (T1Z5), for the low-energy neural- 
network bin (<30keV). (Top panel) Contours of constant 
neural-network output (grayscale solid curves, see text) in the 
plane of the first two principal components. ^®^Cf NRSS and 
NRMS events (blue open circles) and ^^^Ba WBMS events 
(red crosses) are also shown. (Bottom panel) Distributions 
of the “transformed” neural-network output: the distribution 
of NRs from the ^®^Cf data mentioned above (blue solid); 
and the surface-event distribution from the ^®^Ba data men¬ 
tioned above (red dashed). The cut threshold (vertical black 
dashed line) with the WIMP-search data passing (green as¬ 
terisk) and failing (black x’s) the timing cut are also plotted 
along the top of the plot. The transformed neural network 
is used to make the distribution separation more visible; it is 
computed by using the inverse of the neuron response func¬ 
tion: </> : (—00,00) — >■ (0,1) with (j){x) = 1/(1 -I- e“^) 

trix (T{Er', a) was similarly fit using 

<T'^{Er]a) + ( 6 ) 

where the bi{a) are matrices of the free parameters for 
each pair of timing quantities and particle type a. The 
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functional form of the variance was motivated by noting 
that in a simple model of a pulse with a linear rise but 
constant rise time, the probability for a noise fluctuation 
before the pulse rises above the noise is inversely propor¬ 
tional to the slope, and therefore is proportional to the 
inverse of the amplitude (energy). 

The was then formed for every event according to 
the following formula: 

= (7) 

Here ^ is the vector embedding the five timing variables 
for each event and the dependence on and a has been 
left implicit on the right-hand side. 

The surface-event goodness-of-fit variable Xb was con¬ 
structed for each event by the definition 

Xb = min(Xp,Xg), (8) 

where Xg and Xp are the charge-side surface-event and 
phonon-side surface-event goodness-of-fit variables re¬ 
spectively. 

Two restrictions were set in the plane of Xb versus Xn 
that together complete the definition of the 5d-x^surface- 
event rejection cut (see Fig. 9). First the potential WIMP 
events were required to have Xn — where Ci is a 
value for the fth detector, set by requiring that 90% 
of the calibration-data neutrons pass the cut. To dis¬ 
tinguish potential signal events from surface events, it 
was required that Xb ~ Xn — where the index i 

indicates the detector number and the parameter e in¬ 
dicates the event’s energy bin (10-20 keV, 20-30 keV, or 
30-100 keV). 

We parametrize-using adjustable parameters t and b- 
the ionization-yield (y) restriction as 
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fj>nr bdnr ’E V ’E y^nr T t^nr: (9) 

where ynr{Er) and Unr{Er) are the energy-dependent 
mean and standard deviation of the ionization yield for 
NRs, as found using neutron calibration data. The pa¬ 
rameters t and b were required to be the same for all 
detectors and energy bins. The values rji{e), t, and b 
were determined by a simultaneous optimization for all 
detectors and energy bins that maximizes the total SAE 
for eOGeV/c^ WIMPs. 

Optimization of the 5d-x^ timing cut and the NR band 
definition was based on requiring the best overall ex¬ 
pected sensitivity. For a given timing cut (the set {rji}) 
and NR band definition, the expected sensitivity is con¬ 
structed by dividing the 90% Poisson upper limit on the 
expected leakage by the SAE. For an individual detector 
we denote the leakage as Ci{rn) and the SAE as Si{r]i). 
These are smooth functions computed by fitting leakage 
and exposure evaluations at discrete rji using ^^^Ba cali¬ 
bration data. 

Because the cut is defined such that decreasing r]i will 
loosen the restriction, both the leakage and the SAE 
for each detector are monotonically decreasing-leading 


FIG. 9. (Color online) Distributions of calibration data in 
the 5d-x^ parameters for a representative detector (T1Z2). 
(Top panel) ^®^Cf events in the NR ±2cr band (blue open 
circles) are shown with the consistency cut that retains 90% 
of neutrons (vertical black dashed line) in the Xb versus Xn 
plane. (Bottom panel) Events passing the consistency cut- 
NRs from ^®^Cf data (blue open circles) and surface events 
from ^^®Ba data (red crosses)-in the Xb “ Xn versus recoil- 
energy plane. The energy-dependent 5d-x^ timing cut is also 
shown (black solid line) and events above the line pass the cut. 
Note the cut is tighter at low energies and looser at higher 
energies. 

to the condition that the slopes dSildCi for the opti¬ 
mum cut are equal. If not, a unit increase in the leakage 
of a detector with larger slope could be offset by a unit 
decrease in a detector with a smaller slope with a net 
increase in exposure, improving the sensitivity. The tim¬ 
ing cut optimization was done with respect to this slope, 
which parametrizes the {rji} uniquely. 

For the optimum set of timing parameters { 77 ^} the NR 
band definition is selected by choosing the values of t and 
b that optimize the overall expected sensitivity. 

The yield and timing cuts that optimize the expected 
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sensitivity for a 60 GeV/WIMP produced an asym¬ 
metric NR band cut with b = —1.9 and t = 1.8. The 
sensitivity optimization gives a total expected leakage of 
0.5 events and a total SAE of 250 kg days given a WIMP 
mass of 60GeV/c^. The estimated leakage after imblind- 
ing is 1.19 (-1-0.23-0.21)stat events. 


D. Extended analyses 

There has been growing interest in low-mass WIMP 
searches because of some intriguing published results [39, 
40] and the suggestion that the baryon asymmetry is re¬ 
flected in the dark matter sector [32]. While many pre¬ 
vious WIMP searches-guided by the SUSY neutralino 
parameter space-paid much attention to ~100GeV/c^ 
WIMP masses, it is interesting in light of these new re¬ 
sults to examine data with techniques optimized for much 
lower WIMP masses ^lOGeV/c^. For this reason the 
timing-cut constructions considered so far were extended 
with lower thresholds (down to 5 keV for some detectors), 
in order to improve sensitivity to low-mass WIMPs. 

All of the timing cuts presented so far were optimized 
for a 60 GeV jc? WIMP mass, and the corresponding 
analyses were restricted to a recoil-energy threshold of 
10 keV. For events above this threshold the best estimates 
of surface-event background leakage are about one event. 
While the surface-event leakage is worse for the extended- 
threshold analyses, for low WIMP masses sensitivity im¬ 
proves considerably because of the steeply rising WIMP 
spectrum. Simple extensions to the three timing analyses 
were accomplished by lowering the thresholds, confirm¬ 
ing the cut efficiencies below 10 keV, and reevaluating the 
surface-event leakage estimates. The analysis region for 
the extended analyses is approximately 5-15 keV (thresh¬ 
old differs by detector; see below), since recoils of light 
WIMPs (<I0GeV/c^) with energies greater than I5keV 
are very rare. 

For the extended analyses there are several small 
changes to the event selections that help to maximize 
the sensitivity to WIMP masses below lOGeV/c^. 

A careful study of the leakage induced by lowering the 
charge threshold showed that for light WIMPs with sim¬ 
ilar cross sections to the GDMS II silicon result [39] a 
net gain is obtained by lowering the charge threshold to 
3.0 standard deviations above the mean noise value. We 
therefore modified the charge threshold to this value for 
the extended analyses. The recoil-energy threshold was 
then effectively set by the condition that signal events 
have an ionization yield below the mean -3 (t ER band 
line. Of course, signal events were still required to be 
within the mean ±2(T NR band (with the usual slight 
modification for the 5d-x^ analyses). Finally, the signal 
regions were taken to range from threshold up to 15 keV 
since this is the region that is expected to add signifi¬ 
cantly to the low-WIMP-mass exposure. These changes 
are implicit where we use the abbreviation “NRSS” in 
the context of the extended analyses. Most of the timing 


cuts (except the Sd-y^, see below) did not have their deci¬ 
sion boundaries reoptimized for lower WIMP masses; the 
thresholds were simply extended and the changes above 
were made. 

The 5d-y^ method, in addition to lowering the thresh¬ 
olds, was partially reoptimized for a 8GeV/c^ WIMP 
mass. The definitions of the variables (yNAe) 
mained the same. Additionally, the means and covari¬ 
ances were taken to have the same functional depen¬ 
dences given in Eqs. (5) and (6). The optimized yield 
cut was an asymmetric cut, with parameters t = 1.8 
and h = —1.9, the same as the regular 5d-y^ analysis. 
The reoptimized sensitivity for the 8 GeV /<? was not sig¬ 
nificantly better than the standard 5d-y^ analysis, and 
not as good as the sensitivity for the classic extended- 
threshold analysis. Therefore, this reoptimization was 
not carried any further and the 60 GeV j(? optimized ver¬ 
sion is used. 


E. Timing efficiencies and summary 

Energy-dependent efficiencies for the timing analyses 
were computed using neutrons from the ^®^Cf calibra¬ 
tion data. Figure 10 presents the combined efficiencies 
of all cuts (quality, fiducial volume, yield, and timing) 
for each timing analysis. Above a 10 keV threshold, the 
5d-y^ analysis has the best sensitivity and the highest 
SAE for a 60GeV/c^ WIMP (see Fig. 10 caption). The 
classic analysis has the second-best SAE, and provides 
continuity with the original analysis of this data set [33] . 
The neural-network analysis yields the smallest SAE, but 
has the best efficiency just above a 10 keV recoil energy 
threshold. For the extension of the analysis to below 
lOkeV, the classic timing analysis has the best sensitivity 
(though not the highest efficiency at low recoil energies; 
see Fig. 10). As mentioned in Sec. V, the 5d-x^ is our 
“primary” method above lOkeV and the classic is the 
primary extended-threshold method. 

The efficiency function is a necessary ingredient for 
producing the limits and any uncertainty on this function 
is also present in the final limit. The trigger efficiency un¬ 
certainty is ~1% across the energy range; this is mostly 
statistical uncertainty. Our quality cut efficiency is cal¬ 
culated based on baseline noise levels using large event 
populations and so has negligible uncertainty. Efficien¬ 
cies of the other cuts are measured by selecting neutron 
populations in ^^^Gf data and observing the decrease in 
the population by the application of the cuts in bins of 
recoil energy. The results are then fit with an empiri¬ 
cal functional form with a low-energy falloff similar to 
the error function. Fiducial-volume, ionization-yield and 
phonon-timing cuts each contribute about a 5% statisti¬ 
cal uncertainty. This measurement method is, however, 
prone to error due to the fact that neutrons can have 
multiple scatters inside a detector, while WIMPs cannot. 
In the case of the fiducial-volume cuts we used a Monte 
Carlo simulation to find a 5% discrepancy for multiple 
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FIG. 10. (Color online) Total combined efficiencies for all tim¬ 
ing analyses. (Top panel) Efficiencies for the lOkeV threshold 
analyses. The SAEs (60GeV/c^ WIMP) are 219.1 kg days for 
the classic analysis (black solid), 216.4kgdays for the neural- 
network analysis (red dot-dashed), and 262.3 kg days for the 
5d-x^ analysis (blue dashed). (Bottom panel) The efficiencies 
for the extended-threshold analyses with the addition of the 
5d-x^ analysis optimized for a 8 GeV/WIMP mass (ma¬ 
genta dotted). All of the analyses have a total exposure (be¬ 
fore efficiency reductions) of 612.2 kg days. 

scatters, and corrected for it. Overall we assign a conser¬ 
vative 3% systematic uncertainty on the fiducial-volume 
and yield cuts for the effect of multiple scattering. Based 
on these estimates, using the averaged sizes of each effi¬ 
ciency, we expect the total uncertainty on the efficiencies 
to be ^6%. Generally this is negligible on the scale that 
our final limits are presented. A specific study of the effi¬ 
ciencies for the extended-threshold analysis below 6keV 
showed that for the extended-threshold limits this total 
uncertainty becomes ~10% at the lowest WIMP mass 
(6.26GeV/c^) and drops to ^7% by a WIMP mass of 
7GeV/c2. 


VI. UNBLINDING 

A blinding technique was used to avoid bias in the set¬ 
ting of data selection cuts. The current work used what 
can be referred to as “hidden signal box” analysis [68] 
which is common in rare-event searches where the signal 
region is known a priori. The same data were analyzed 
previously [33], but since all the data were reprocessed 
with an upgraded charge reconstruction algorithm, and 
the cuts were optimized solely based on calibration data 
and distributions outside the newly masked signal region, 
this is a good approximation of a hidden signal box anal¬ 
ysis. Signal events were hidden by removing the single¬ 
scatter events in the NR yield band. This technique is 
effective for removing bias in the timing-cut preparation, 
but restricts the information that can be used for estimat¬ 
ing surface-event leakage before unblinding has occurred. 
For this reason the leakage estimates were done before 
and after the unblinding, but consistency between the 
two methods is checked. Including the post-unblinding 
version using the NR singles that fail the timing cut made 
the final estimate more robust. 

Upon unblinding, the following number of events pass 
all cuts above lOkeV: zero for the 5d-x^ timing analy¬ 
sis; two for the classic timing analysis; and one for the 
neural-network analysis. One of the two candidates from 
the original analysis is a candidate in both the classic 
timing and the neural-network analyses [33] . The second 
candidate from that analysis, whose poor charge-pulse 
fitting prompted this reanalysis (see Sec. Ill), failed all 
timing analyses by a substantial margin. Information 
about the two passing candidate events is shown in Ta¬ 
ble II. Figure II shows the location of these events within 
the signal region of the classic analysis. 

Detector Recoil energy [keV] Yield Analysis 
T1Z5 12)30 0.33 NN, Glassic, 2010 

T2Z3 10.81 0.33 Glassic 

T3Z4 15.35 0.26 2010 


TABLE II. Information about the three WIMP candidate 
events above lOkeV. Under “Analysis,” “NN” refers to the 
neural-network timing analysis and “2010” refers to the origi¬ 
nal analysis of these data [33]. Note that no candidate events 
are observed for the 5d-x^ analysis. 

Our extended-threshold analyses gave a wide range in 
terms of the number of candidate events. The classic 
analysis had six candidates and the neural-network anal¬ 
ysis had 16. No events were observed in the extended 
version of the 5d-x^ analysis. A description of the can¬ 
didate events for our primary extended-threshold anal¬ 
ysis (the “classic” cut) is given in Table III. The large 
number of candidate events in the neural-network anal¬ 
ysis is attributed to an increased leakage of anomalously 
low-ionization events, those that are normally below the 
ionization threshold in the 10 keV analysis [62] . This in¬ 
creased leakage of essentially zero-ionization events was 
not expected and is presumably due to a bias in the train- 
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timization to exposure at a 60GeV/c^ WIMP mass. 


Detector Recoil energies [keV] 

T1Z5 

3.45, 5.73, 12.30 

T2Z3 

10.81 

T4Z4 

7.56 

T4Z5 

7.25 


TABLE III. Information about the six WIMP candidate 
events for the “classic” extended-threshold analysis. 

The post-unblinding leakage estimates can be found in 
Table IV along with the analysis energy ranges, candi¬ 
date numbers, exposures and WIMP-mass optimization 
assumptions. 


VII. ESTIMATED BACKGROUNDS 
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FIG. 11. (Color online) Distributions of WIMP-search events 
in the detectors containing the two candidate events for the 
classic timing analysis: T1Z5 (top panel) and T2Z3 (bottom 
panel). The normalized yield is the distance from the average 
NR yield in units of the standard deviation of the yield distri¬ 
bution; the black horizontal lines on the right side of the plot 
indicate the ±2 ct NR band. The normalized timing parameter 
is the standard timing parameter (tdei+r) minus the value of 
the cut boundary. WIMP-search events (red crosses) to the 
left of the timing cut (vertical black dashed line) pass all cuts 
except phonon timing and yield; events to the right of that 
line pass the timing cut, and those that are also in the NR 
band (solid black lines) pass the yield cut. The highlighted 
event (blue open circle surrounding red cross) in the top panel 
is candidate 1 and the highlighted event in the bottom panel 
is candidate 2 (see Table II). All the events shown are in the 
recoil energy range 10-100 keV. 

ing set of the neural network. It has negligible contri¬ 
bution above the original 4.5 standard deviation charge 
threshold. The low number of events in the 5d-x^ analy¬ 
sis is not unexpected because that cut has an independent 
energy bin at 10-20 keV, where the cut is rather stringent 
because of increasing leakage at low energy and the op¬ 


A. Surface electron-recoil background 

Surface electron recoils originate from several sources: 
1) particles emitted from /3 emitters contaminating the 
surfaces of the detector and the material around it (no¬ 
tably ^^°Pb), 2) photo-electrons emitted from material 
neighboring the detector through Compton scattering or 
the photoelectric effect, and 3) photons that interact in 
the detector within a few microns of the surface. Pho¬ 
tons from category 3) can be low-energy xrays or high- 
energy photons that Compton scatter in the detector. 
Past studies showed that the dominant contributions are 
from ^^'^Pb and photon-induced backgrounds, which con¬ 
tribute approximately equally. No other sources were 
found to be statistically significant. 

The expected number of surface events leaking into 
the WIMP signal region was calculated using the num¬ 
ber of single scatters in the NR band that are rejected by 
the timing cut and the surface-event rejection efficiency. 
The surface-event rejection efficiency was estimated from 
three independent event sets and combined to improve 
accuracy. To reduce systematic uncertainties, leakage es¬ 
timates were calculated on a detector-by-detector basis 
(index i below) and where possible the relevant event set 
was separated into bins (index j below) of energy and 
approximate event position (see “face” bins below). In 
every case the leakage can be expressed as: 

= ( 10 ) 

where n is the total expected number of surface events 
leaking into the WIMP signal region. The symbol W 
is the number of NRSS events in the NR band rejected 
by the timing cut for the zth detector. and rriij are 
the number of multiples in (or around; see below) the 
WIMP signal region failing and passing the timing cut, 
respectively. The Sij are the fractions of Ni in subset j, 
which are calculated using the surface-event multiples in 
the NR band for the WIMP-search data. 
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Method 

Energy range [keV] Exp. leakage Candidates WIMP 

mass [GeV/c^j 

SAE [kg days] 

Classic 

10-100 


2 

60 

220 

Neural Network 

10-100 

0.87lE]-Jt 

1 

60 

216 

5d-x" 

10-100 

1.19l°J? 

0 

60 

250 

Classic ext. 

5-15 


6 

60 

186 

Neural Network ext. 

5-15 


16 

60 

190 

extended 5d-x^ 

5-15 

> 0.97t[;;« 

0 

60 

202 

extended 5d-x^ 

5-15 


0 

8 

3.78“^ 


SAE depends on an assumed WIMP mass [see Eq. (4)]; we use the optimization mass in all cases. For this reason SAE is only 
comparable in situations of common optimization mass and signal-region energy range 

TABLE IV. Expected leakage and exposure statistics for all of the surface-event rejection methods described in this work. For 
the extended-threshold analyses, we quote 5keV as an approximate lower limit on the signal region. The actual threshold 
depends on the detector and is set by the crossing of the ER band limits and the charge threshold curves (see Sec. VD). The 
symbol > is used to indicate lower limits on the expected leakage. In those situations the event sets that are typically used to 
estimate the leakage do not have events all the way down to the signal-region threshold (see text). 


The 14 detectors used in this analysis were split into 
two detector sets according to their positions in the 
tower: the 12 interior detectors, and two “endcaps,” on 
the top or bottom of a stack of six detectors. Equa¬ 
tion 10 was applied to each detector set separately be¬ 
cause endcaps require additional systematic uncertainty 
corrections (discussed below). The and rriij are de¬ 
termined independently for three different event sets: 1) 
NRMS in WlMP-search data, 2) WBMS in WlMP-search 
data, and 3) WBMS in ^^^Ba calibration data. For rep¬ 
resentative surface-event populations drawn from outside 
the NR band, the index j specifies one of six subsets con¬ 
structed from two detector faces (see Sec. IV B on face 
tagging) and three energy bins (10-20 keV, 20-30 keV, 
and 30-100 keV). In the case where and rriij were 
drawn from the NRMS WlMP-search data, the six event 
subsets were combined for each detector because of the 
low statistics, reducing the index set for j to one ele¬ 
ment. Each of the three independent estimations of rriij 
and Mij are used with Eq. (10) to obtain three estimators 
for n, which were then statistically combined to produce 
the final leakage estimate. 

With these estimators, Monte Carlo (MC) simulations 
and Bayesian inference were used to estimate the surface- 
event background and its uncertainties after unblinding 
the data [66]. Instead of simulating the posterior distri¬ 
bution of n (Eq. (10)), for each event set, MC simulations 
are run for the posteriors of the individual Poisson counts 
rriij and Mij, and the multinomial fractions Sij. Jaynes 
priors, p{p) oc with c « — 1, were chosen for the Pois¬ 
son distribution. This prior is generally considered an 
“objective prior,” having the advantage of ensuring in¬ 
variance of statistical inference under transformations of 
the Poisson mean [69]. A uniform prior was used for the 
multinomial distribution. The posterior of n was then 
calculated using Eq. (10) with the simulated posteriors 
of each component. 

Systematic uncertainties of different origins were es¬ 
timated independently, added in quadrature, and then 
incorporated into the posterior of the leakage with a pro¬ 


file of the standard normal distribution. Systematic un¬ 
certainties from two sources were estimated, for all the 
detector and data-set combinations, including the choice 
of the prior exponent c and the difference between singles 
and multiples. 

End-cap detectors can have their multiples tagged only 
on one side. This biases the counts of singles upward in 
these detectors. Together with the low number of sur¬ 
face events in the data, it is challenging to estimate the 
surface-event background. As a result, a conservative 
estimate where we allow the Sij to be biased higher for 
endcaps-the detectors with the worst leakage-is used. As 
a double check, the component leakages for each detec¬ 
tor obtained by the Bayesian approach were compared to 
those obtained by the frequentist approach [70]. There 
is good agreement between the two. 

The final surface-event leakage estimates incorporat¬ 
ing both systematic and statistical uncertainties are dis¬ 
played in Table IV for all timing cuts presented in this 
work. The analyses with the best expected sensitivities 
give leakage estimates of 1.19 (-1-0.23-0.20) events (5d- 
for the lOkeV threshold analysis) and > 1.48 (-1-0.20- 
0.20) events (classic extended threshold). The leakage 
estimate for the extended-threshold analysis is a lower 
limit because our standard background sample does not 
extend all the way down to the signal-region threshold. 
WlMP-search multiple scatters below the ER mean -5 ct 
line and above the NR mean -|-2cr line (WBMS) is one 
of our typical background samples, but because the ER 
mean -ba and NR band lines cross there is an implicit 
energy threshold in this background sample. The signal 
region, however, is taken to be the NR band above the 3cr 
charge threshold and below the ER mean -3cr line. This 
leaves a small region (see Fig. 5) that is unaccounted for 
in the background estimates. The estimate is neverthe¬ 
less useful because 1) the implicit background-estimation 
threshold can be as low as ^5 keV, very near the signal- 
region threshold in most cases, and 2) the limits are cal¬ 
culated without background subtraction so a conspicuous 
rise in the WIMP upper limit would exist if there were a 
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large leakage present in the unaccounted region. 


B. Cosmogenic background estimate 

To estimate the cosmogenic neutron background, a 
combination of simulation and data was used [71] . Parent 
muons that intersect the CDMS II 5-cm-thick scintilla¬ 
tor panels are vetoed at nearly 100% efficiency. Our data 
contains only a small number of neutrons that accompany 
identified muons, so a high-statistics simulation was used 
to estimate how many veto-coincident neutrons produce 
NRs in the Ge ZIP detectors, compared with those that 
are not accompanied by veto activity. 

Simulation muons generated using the MUSUN [72] 
program based on slant-path data from Soudan2 an¬ 
gular muon flux measurements [73] were used as input 
to a GEANT4 MG model of the Soudan GDMS II ex¬ 
perimental setup. These muons have a mean energy 
of ^215 GeV and azimuthal and zenith distributions 
characteristic of the overburden of the Soudan site and 
depth (2090m.w.e.). The simulated muons were gen¬ 
erated with the appropriate angular distributions and 
spectra on a five-sided parallelepiped (no floor). They 
were then propagated via GEANT4 through 10 m of rock 
into the Soudan hall. All secondaries and the parent 
muon(s) were tracked through a complete geometry of 
the GDMS II shielding and detector towers. In the sim¬ 
ulations, the GDMS setup is located asymmetrically 2 m 
from one wall of the 8 m east-west cavern dimension to 
reproduce albedo effects. The statistics correspond to 66 
live years. A multiple is defined in the simulation output 
as an NR in a Ge detector accompanied by energy depo¬ 
sition of any sort above 2 keV in any other ZIP detector. 

In an effort to balance the statistical and systematic 
uncertainties, three complementary estimators for the to¬ 
tal number of unvetoed single-scatter events were con¬ 
structed. For the first two estimators, the simulation 
was used to calculate a veto ratio (unvetoed to veto- 
coincident NRs) that was then normalized to the veto- 
coincident data from the WIMP-search measurements to 
establish the background estimate. The second estimator 
makes use of the higher-statistics multiple-scatter sam¬ 
ple, whereas the first uses only single scatters. The third 
estimator uses only simulated unvetoed single-scatter 
events, scaled to the correct experimental live time, and 
accounts for the detector efficiencies for each timing anal¬ 
ysis. 

The veto ratio is 0.008 ±0.003 for single scatters and 
0.003 ± 0.001 for multiple-scatter NR events in the sim¬ 
ulation. One might expect a smaller veto ratio for 
multiple-scatter events if they are taken to represent 
more pervasive showers; excess energy from those showers 
will often be detected by the veto system even when the 
parent muon misses the scintillator panels. The weighted 
average for all events gives a veto ratio of 0.004 ±0.001. 

There are 14 Ge NR events in our veto-coincident data, 
of which three are single scatters. This is in agreement 


with the MC which predicts that ~26% of the vetoed 
events should be singles, and the rest multiples. Of the 
vetoed NR events, five pass the 5d-x^ timing cut and 
four pass the classic or neural-network timing cuts. For 
the data-driven estimators, we normalize to those passing 
the timing cut, rather than introduce additional system¬ 
atic uncertainty by applying an average timing-cut effi¬ 
ciency to all 14 detectors in the simulation. However, the 
numbers obtained either way are consistent, with trade¬ 
offs between systematic uncertainty and statistical uncer¬ 
tainty in each case. We decided to quote the simulation- 
driven estimate because it is consistent with the others 
and offers the best statistical uncertainty. The system¬ 
atic uncertainty is taken as the spread between all the 
estimators. 

Based on the 16 unvetoed NR singles over 66 live 
years of simulation data, the 5d-x^ timing-cut cosmo¬ 
genic neutron background is 0.021 ± 0.008stat ± O.OOQgys 
events for the lOkeV analysis. The classic and neural- 
network cuts give 0.019 ±0.007stat±0.01sys events and 
0.018 ±0.007stat±0.01sys events respectively for the 
lOkeV threshold case. 

The estimates were also made for the extended- 
threshold analyses by using the energy range 2-20 keV 
to approximate the extended-threshold energy ranges. 
Based on 12 unvetoed NR singles in that range over 
66 live years the 5d-x^ extended timing-cut cosmogenic 
neutron background is 0.009 ± O.OOdgtat ± O.OOlgyg events, 
the classic is 0.012 ± O.OObgtat ± O.OOdgyg events, and the 
neural network is 0.014 ± 0.006stat ± O.OOOgys events. 


C. Radiogenic background estimate 

The radiogenic neutron background was also estimated 
using a GEANT4 simulation of the CDMS II tower con¬ 
figuration. Neutrons originating from the Th and U de¬ 
cay chains were simulated for each material in the setup 
(lead, polyethylene, copper). The primary energy spec¬ 
tra of the neutrons was generated using the SOURGES4G 
package [74-76], which computes the neutron spectra due 
to spontaneous fission and (a, n) reactions of alphas from 
the full decay chains within the matrix material. The 
neutrons were propagated through all materials, eventu¬ 
ally creating NRs in the Ge and Si detectors. The single- 
and multiple-scatter rates were tabulated for each detec¬ 
tor to create high-statistics files labeled by their contami¬ 
nant source. The files were then weighted by the contam¬ 
ination level (see below) and normalized to the amount of 
material present to determine the final background rate. 
The source with the highest single-scatter contribution 
is 238 u copper cryostat enclosures (cans), followed 

closely by U and Th contaminants in the lead shielding. 

The contamination levels were determined by a sep¬ 
arate 7 simulation, again using GEANT4 and the same 
CDMS II geometry. In this simulation, gammas from the 
^^^Th and ^38^ decay chains, and ®°Co were gen¬ 
erated from inside the shielding components and tower 
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structures, and from radon daughters on surfaces. In or¬ 
der to reproduce the fiducial-volume cut already applied 
to the data, the location of the energy deposition and the 
electric field map in the detectors were used to deduce the 
fraction of charge collected on the inner and outer elec¬ 
trodes. The event was cut if it produced a measurable 
signal on the outer electrode. A charge threshold was 
implemented in the simulation data by using the exper¬ 
imental inner-electrode charge thresholds (in keVee) and 
applying them to the inner-electrode energy derived from 
the above procedure. 

To obtain the contamination estimate a minimiza¬ 
tion was performed to fit the 43 sources considered (see 
Fig. 12 for a schematic of the setup and Table V for a 
breakdown of source locations) to the ER data above 
15 keVee, avoiding the ~10keVeeGe activation peaks. 
Some representative values for this contamination are 
4 mBq/kg for both ^^^Th and in the inner polyethy¬ 
lene shield and 6-7 mBq/kg for ^^^Th and in the 
tower 5 copper, the tower with the highest contamina¬ 
tion. The comparison between the ER spectrum and the 
sum of MC results using the best-fit contamination levels 
is shown in Fig. 13. 


CDMS Shielding and Veto 



FIG. 12. (Color online) Schematic of the CDMS II shield¬ 
ing configuration. The outermost vacuum can surface is 
shown near the center and radially outward the layers are: 
inner polyethylene (green), inner lead (light grey), outer lead 
(darker grey), outer polyethylene (green), and muon-veto 
scintillator counters (blue) with attached light guides and 
phototubes (white and black). Cooling of the detectors is 
achieved through the “cold stem” from a dilution refrigera¬ 
tor (dark blue), while cabling passes through the “electronics 
stem” on the other side of the setup. 

Because of nearly degenerate fits to the spectrum, 
there is uncertainty in the final contamination values. 
The degree of uncertainty was estimated by bracket¬ 
ing the largest changes that produce acceptable fits in 
the gamma MC and translating these into the resulting 
changes in the neutron rate. 

In the 14 Ge detectors used for the WIMP-search anal¬ 


ysis, the expected raw radiogenic neutron single-scatter 
background rate was found to be (1.15 ± 0.14) x 10“^ 
events/kg/day in the 10-100 keV energy range. The 
statistical uncertainty is negligible because of the large 
number of events simulated while the systematic uncer¬ 
tainty listed is due to the uncertainty in the contami¬ 
nation values and locations. This singles spectrum was 
then convolved with the final neutron cut efficiencies from 
the three timing analyses (see Fig. 10). The radiogenic 
spectrum-averaged efficiencies are: 35.3% for the clas¬ 
sic, 34.9% for the neural net, and 40.4% for the 5d-x^ 
analysis. For the raw exposure of 612.2 kg days, the final 
radiogenic neutron event backgrounds for the three anal¬ 
yses are: 0.025 ± O.OOOlstat ± 0.003sys events for the clas¬ 
sic and neural network, and 0.028 ± O.OOOlgtat ± O.OOdgyg 
events for the 5d-x^ analysis. 

For the extended analyses the estimates were made 
with the approximate energy range 2-20 keV. The esti¬ 
mates were: 0.0131 ± O.OOOlstat ± 0.0019sys for the clas¬ 
sic, 0.0148 ± O.OOOlstat ±0.0021sys for the neural net¬ 
work, and 0.0105 ± O.OOOlstat ± 0.0015sys for the 5d-x^. 

D. background estimate 

^^°Po decays via the a-decay process. This leaves an a 
particle of 5.4 MeV and a recoiling (<100 keV) ^°®Pb nu¬ 
cleus in the final state. ^^°Po is in the chain 

and is part of possible post-^^^Rn plate-out contamina¬ 
tion. It is also the first post-^^^Rn a-decayer after the 
long-lived ^^°Pb. These facts make this decay likely to 
be an important long-lived contamination, and because 
a recoiling ^°®Pb nucleus could have produced an energy 
deposition below ~100keV with an ionization yield con¬ 
sistent with being an NR, it makes sense to evaluate this 
background very seriously. 

Because of the GDMS II detector geometry (see Fig. 1) 
much of the surrounding surface of an interior detector 
is comprised of another detector. ^^°Po decay events in 
which the decay occurs in one detector and the ^°®Pb re¬ 
coil is registered in an adjacent detector cannot contam¬ 
inate the signal region of the latter detector because of 
the clear 5.4 MeV a deposition in the former. Therefore 
the most important component of the background comes 
from surfaces that are uninstrumented or those that are 
adjacent to uninstrumented surfaces. Another important 
point is that since the alphas from this decay give such 
a clear signature, and the decay is a two-body decay, 
an obvious way to estimate the number of unaccompa¬ 
nied (single) ^°®Pb events is to estimate the number of 
unaccompanied a events from the decay. The angular 
distribution is isotropic and any surface that can observe 
a single a is approximately equally likely to observe a 
single ^°®Pb recoil. This estimation of the number of 
single ^°®Pb events was carried out, and using rough es¬ 
timations for the passage fractions of these events when 
subjected to the other analysis cuts, we expect approx¬ 
imately 0.187 signal-region ^°®Pb events over the whole 
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Contamination 

Outer Lead 

Inner Lead 

Inner Polyethylene 

Cans 

Innermost Can Surface 

Outermost Can Surface 

Towers 


/ 

- 

- 

/ 

/ 

/ 

/ 

40k 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

232rpj^ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

238-q 

/ 

/ 

/ 

/ 

/ 

/ 

/ 


- 

- 

- 

- 

- 

/ 

- 


TABLE V. Sources used for the radioactive contamination fitting procedure. The /indicates that the respective source was used 
in the fitting procedure for the location given. The cans are the nested copper cold stages of the cryostat, where surface sources 
were simulated on both the innermost and outermost surfaces. There are five separate sources for the copper components of 
the five detector towers. No simulated sources were placed in the outer polyethylene or the veto panels because they are behind 
too much shielding to yield a measurable amount of background. 



FIG. 13. (Color online) Measured event-summed energy across all detectors (black solid line) compared to simulations of 
contaminants from common locations: lead shielding (blue diamonds), inner polyethylene shielding (magenta crosses), copper 
cans and tower components (cyan filled circles), and outermost can surfaces (brown stars). The sum of all simulated components 
(red inverted triangles) is also plotted and is in good agreement with the measured data. 


exposure. 

Since the recoil estimates were inferred from a 

counting, and the recoils can come at different 

energies depending on how deeply the ^^*^Po parent is 
embedded into the originating surface, we do not have 
a good specification of the energy distribution of such 
events. Therefore, to be conservative we can use the same 
estimate for the lOkeV and extended-threshold analy¬ 
ses. We expect 0.187± O.OlSstat ± O.lSTgys events over 
the whole exposure; and have assigned a 100% system¬ 
atic uncertainty to account for the roughness of this esti¬ 
mate. This background estimate is clearly subdominant 
with respect to the surface-event background estimates. 


but is larger than the cosmogenic or radiogenic neutron 
background estimates [62]. 


VIII. RESULTS 

The background estimates for the primary 10 keV and 
extended-threshold analyses are summarized in Table VI. 
While the background estimates can be used to interpret 
the overall results of the experiments they do not directly 
modify the limit curves (see below). 

Results of direct WIMP-search experiments are usu¬ 
ally summarized as upper limits or signal contours in the 
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Background 5d-x‘^ (lOkeV) Classic (extended) 

Leakage 

1.19 ±0.22 

>1.48±0.20 

Cosmogenic 

0.021 ±0.012 

0.012 ±0.006 

Radiogenic 

0.028 ±0.004 

0.013 ±0.002 

206pb 

0.19±0.19 

0.19±0.19 

Total 

1.43 ±0.30 

>1.69±0.28 


TABLE VI. A summary of the background estimates from 
the primary analysis methods. Statistical and systematic un¬ 
certainties have been added in quadrature and to compute 
the total backgrounds all component uncertainties were also 
added in quadrature. 

WIMP-nucleon cross section versus WIMP-mass plane. 
Yellin’s optimum interval method [77] allows derivation 
of an upper limit on a signal rate in cases with unknown 
background. While the backgrounds in our signal region 
are not completely unknown, this is a conservative ap¬ 
proach to setting upper limits on a possible signal. This 
presentation also requires assumptions about the WIMP 
distribution in the galactic halo, the type of interaction 
between WIMPs and nucleons, and the nuclear form fac¬ 
tor for the interaction. The velocity distribution was as¬ 
sumed to be Maxwellian and was parametrized by the 
rotational velocity at infinite radius and corrected for the 
finite galactic escape velocity [49], which is taken to be 
544 km/s [78]. A WIMP mass density of 0.3 GeV/c^/cm^ 
was used for historical reasons, making the computed lim¬ 
its comparable to similar publications [49]. Some recent 
astrophysical measurements indicate different values [79] ; 
a correction for such deviations is a simple multiplica¬ 
tive factor and can be easily applied by the reader [80]. 
A most probable WIMP velocity of 220 km/s was used 
along with a mean circular velocity of the Earth with 
respect to the Galactic center of 232 km/s. The WIMP 
interactions were assumed to be spin independent and the 
Helm form factor was used [49] for a natural Ge isotopic 
distribution. 

The comparison of the present 10 keV 5d-x^ and clas¬ 
sic extended results with other published limits and sig¬ 
nal contours is shown in Fig. 14. In the figure, our 
10 keV 5d-x^ limit is combined with the CDMS II five- 
tower exposure acquired before July 2007, resulting in 
a limit that summarizes the full (and final) GDMS II 
high-threshold sensitivity. The CDMS II/EDELWEISS 
combined limit [35] is also shown for comparison. Above 
~I00 GeV/c^ WIMP mass the combined limit is compa¬ 
rable to our CDMS II combined result owing to the good 
efficiency-averaged exposures of both of the experiments 
in the relevant energy ranges. 


A. Limit cross-checks 

To gain insight into the effect of timing-cut figures of 
merit (efficiency, SAE, leakage) on WIMP-search results, 
we have constructed final limits for all of the timing-cut 
constructions in this work. Although the SAE and ex¬ 
pected leakage were exclusively used to choose the pri- 



FIG. 14. This figure compares the main results from 
this analysis (the current Sd-y^ analysis combined with prior 
CDMS II exposures, black solid and the classic extended- 
threshold analysis, black dashed) with previously published 
results (all limits are at 90% C.L.): Darkside-50 [81] (orange 
triple-dot-dashed) XENONIOO [82] (blue triple-dot-dashed); 
LUX [83] (red dot-dashed); SuperCDMS low-threshold [84] 
(green dotted); CDMS II/EDELWEISS combined [35] (pur¬ 
ple long dashed); CRESST-II [85] (magenta long-dot-dashed); 
CDMS II Si [39] (90% and 68% C.L. contours, blue, with the 
best-fit point marked with a black dot); DAMA/LIBRA [43, 
86] (3(7 region, light red) and CoCeNT [40] (90% C.L., 
brown). 

mary timing cuts, comparing all of the limits in this way 
gives cross-checks on how other parameters (efficiency, 
signal-region events) affect the reach of the experiment. 

The top part of Fig. 15 shows the limits derived using 
the optimum interval method for the different timing cuts 
with 10 keV thresholds. Relative to the original publica¬ 
tion, the ionization-based fiducial-volume and phonon¬ 
timing cuts have improved efficiencies in the analysis 
reported here, leading to the improved exposure and 
more stringent limits. In terms of the overall spectrum- 
averaged detection efficiency for a 60 GeV/WIMP, the 
classic timing-cut strategy (see Sec. V A) shows a 12% im¬ 
provement over the previously published version, about 
half of which can be attributed to a reoptimization of the 
ionization-based fiducial volume following the data repro¬ 
cessing. A similar improvement was seen for the neural- 
network timing-cut analysis (described in Sec. VB). The 
largest improvement of 29% in overall SAE efficiency was 
achieved for the Sd-y^ analysis (described in Sec. VC), 
owing primarily to an increased timing-cut efficiency in 
the 15-90 keV energy range (see Sec. VE). 

For the 10 keV-threshold analyses, the 5d-y^ sets the 
most stringent limit at a 60 GeV j(? WIMP mass, while 
the neural-network timing cut results in stronger limits 
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at and below 10 GeV/c^, an important region for further 
study [39]. The 5d-x^ set weaker limits for low-mass 
WIMPs because the cut was set to maximize sensitivity 
to a WIMP with mass 60 GeV/c^. This combined with 
the fact that the 5d-x^ method could set a very tight cut 
at low energies (it used an independent 10-20 keV bin 
whereas the other analyses were less granular) produced 
a poorer WIMP efficiency toward low recoil energies de¬ 
spite the lower expected background leakage at those en¬ 
ergies, but excellent sensitivity for high WIMP masses. 

The lower part of Fig. 15 shows the extended limits 
in the low-WIMP-mass region. Each extended analysis 
constrains the 8-10 GeV/c^ mass region more strongly 
than the higher-threshold analyses, and the classic tim¬ 
ing cut produces the strongest limit near the silicon- 
detector analysis best-fit point of Mw =8.6 GeV/c^ and 
asi =1.9x10“'^^ cm^ [39]. The extended limits are also 
compared to the previous low-threshold CDMS II results, 
which did not use a timing cut [34] . That analysis has a 
larger exposure toward lower recoil energies which ac¬ 
counts for the stronger limit set below a ^7 GeV/c^ 
WIMP mass. The classic analysis presented here has 
a stronger limit by a factor of approximately 2.7 at a 
WIMP mass of ~8.6 GeV/c^. 


IX. CONCLUSION 

The reprocessed data did not produce significant 
changes in the number of signal-region events, indicat¬ 
ing that uncertainties applied in the original processing 
of the CDMS II data set [33] were robust. All three sets 
of higher-threshold timing cuts produced similar limits, 
with small differences consistent with their correspond¬ 
ing exposure-optimization procedures. For example, the 
5d-x^ analysis has a high efficiency at moderate recoil 
energies (30-60 keV), but has a stringent timing cut at 
lower energies. It is well suited to provide the strongest 
limits at high WIMP mass (>60 GeV/c^), but will pro¬ 
duce fewer low-energy signal-region events. On the other 
hand, the classic analysis at 10 keV threshold shows a 
slight weakening of the 90% C.L. limit for WIMP masses 
below about 18 GeV/c^, where sharp increases in the 
limit curves indicate systematics near threshold. The 
neural-network timing cut has been identified as a ro¬ 
bust method with the highest signal efficiency at low en¬ 
ergies (see Fig. 10) and good sensitivity at lower WIMP 
mass (< 9 GeV/c^). The classic extended-threshold limit 
rules out about half of the silicon 68% C.L. region ob¬ 
tained in a previous CDMS II publication [39]. This 
indicates that the low-threshold results from the Ge de¬ 
tectors are marginally compatible with the Si-detector 
measurements taken during the same data period, under 
standard assumptions. Comparisons of such direct de¬ 
tection results on different nuclei will be a powerful tool 
for understanding WIMP dark matter both in terms of 
the fundamental WIMP interactions [87, 88] and possible 




WIMP mass [GeV/c^] 

FIG. 15. (Color online) Experimental upper limits (90% confi¬ 
dence level) derived from each of the analyses presented in this 
work compared with the originally published [53] (black dot¬ 
ted) limits. The CDMS II Si contour is shown with the best-fit 
point marked with a black dot (WIMP mass of 8.6 GeV/c^ 
and WIMP-nucleon cross section of 1.9x10“'^^ cm^) [39]. 
(Top panel) The 10 keV threshold analyses. The Sd-y^ limit 
(blue dashed) is the “primary” high-threshold result to be 
quoted from this work. The neural-network and classic limits 
are shown as red dot-dashed and black solid lines respectively. 
(Bottom panel) The extended-threshold limits, focused on the 
lower WIMP-mass region. The same color code applies ex¬ 
cept that all of the analyses from this work correspond to the 
extended-threshold versions. The classic limit (black solid) 
is the “primary” extended-threshold result to be quoted from 
this work. The extended 5d-y^ limit shown corresponds to the 
timing-cut optimization assuming a 60GeV/c^ WIMP mass. 
For comparison, the previous CDMS II low-threshold limit is 
shown [34] (green triple-dot-dashed). 

backgrounds. 
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